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Abstract 

We consider the three-nucleon system within an effective theory with contact interactions at 
leading order in the large scattering length. We calculate the charge form factor of the triton at 
low momentum transfer and extract the triton charge radius. At this order, no two-body currents 
contribute and the calculation can be performed in the impulse approximation. We also comment 
on the power counting for higher orders. The requirement of a three-body force for renormalization 
of the three-nucleon system explains the previously observed correlation between the triton binding 
energy and charge radius for different model potentials. 
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I. INTRODUCTION 



There has been a considerable interest in physical systems with large scattering length 
recently. This interest was stimulated by the effective field theory program in nuclear physics 
0,0] and by the experimental realization of Feshbach resonances with trapped atoms, leading 
to tunable interactions in atomic systems 0, 0, HI • The scattering of two particles with short- 
range interactions at sufficiently low energy is determined by their S-wave scattering length 
a. If a is much larger than the typical low-energy length scale £, which is given by the 
range of the interaction, the system shows universal properties. The simplest example is 
the existence of a shallow two-body bound state if a is large and positive, but there are 
many more, including the effects of a limit cycle 0, 0] and the Efimov effect |H 0] in the 
three-body system. For a recent review of the universal properties of iV-body systems with 
short-range interactions up to iV = 4, see Ref. (lp| . 

The best known example of a nuclear system with a large scattering length is the two- 
nucleon (NN) system. There are two independent S-wave scattering lengths that govern the 
low-energy scattering of nucleons. The scattering lengths a s = —23.5 fm and a t = 5.42 fm 
describe iViV scattering in the spin-singlet ( 1 S'o) and spin-triplet ( 3 Si) channels, respectively. 
Both scattering lengths are significantly larger than the natural low-energy length scale 
I ~ l/m T « 1.4 fm, while the effective ranges are of the same order as 1} As a consequence, 
the description of few-nucleon systems in an expansion in l/\a\ is useful. It has successfully 
been applied to various two-, three-, and four-nucleon observables (See Refs. [H EH 0, E3| 
and references therein). 

The most well-known example of a universal feature of the three-nucleon system is the 
Phillips line If the predictions of different nucleon-nucleon potentials for the triton 

binding energy B t and the spin-doublet neutron-deuteron scattering length a^ 2 ^ are plotted 

against each other, they fall close to a line. This correlation between B t and a^d^ can not 
be understood in conventional potential models. However, it immediately follows from 
universality if the large iViV scattering lengths are exploited within an expansion in £/\a\ 
[Til 1^ | . If corrections of order £/\a\ are neglected, all low-energy 3-nucleon observables 
depend only on the spin-singlet and spin-triplet scattering lengths a s and a t and the three- 
body parameter L 3 . Since the NN potentials reproduce the scattering phase shifts, they 
all have the same scattering lengths. However, the off-shell behavior of the potentials is not 
constrained by the phase shifts and in general is different for each potential. This difference 
is captured in the three-body parameter L3. The various potential model calculations must 
therefore fall close to a line which is parametrized by the parameter L3. A similar universal 
feature of the four-nucleon system is the Tjon line: an approximately linear correlation 
between the triton binding energy B t and the binding energy of the a-particle B a . This 
correlation was discovered by Tjon [17[ using simple separable interactions, but also holds 
for modern phenomenological potentials The origin of this correlation was explained in 
Ref. |13j] from the absence of a four-body force at leading order in the pionless effective field 
theory. As a consequence, the Tjon line is also parametrized by the parameter L 3 . However, 
the universal properties resulting from the large scattering length are not restricted to purely 
hadronic observables. For example, a correlation between the triton binding energy and 
charge radius was observed by Friar, Gibson, Chen, and Payne 1 1 ( .j . 



In the following, we will simply use a as a generic symbol for a s and a t if no distinction is required. 
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This work is another step towards the treatment of external currents in the three-nucleon 
system within the pionless effective theory. We calculate the charge form factor of the 
triton at leading order in £/\a\ and extract the charge radius. Furthermore, we study the 
correlation between the triton binding energy and charge radius and relate it to the three- 
body parameter L 3 that parametrizes the Phillips and Tjon lines. For a recent calculation 
of the cross section for the reaction nd — > 3 H7 at energies below 200 keV see Ref. jiol ]. 



II. EFFECTIVE THEORY 



Effective theories are a powerful tool to calculate low-energy observables in a systematic 
fashion. They are ideally suited to exploit a separation of scales such as the one between \a\ 
and I. At sufficiently low-energies, one can use an effective theory with contact interactions 
only. In this paper, we work in a quantum-mechanical framework and construct an effective 
low-energy interaction potential. In the following, we briefly review this construction. For 
a more detailed discussion see Ref. |2l| . 

The bare effective potential generated by short-range contact interactions can be written 
down in a momentum expansion. In the S-wave sector of the two-nucleon system, it takes 
the general form 



(k'\V bare \k) =V S 



2.2 



(k 2 + k l2 )/2 



Vt 



2.2 



(k + k u )/2 



+ ... 



where the dots indicate momentum dependent terms that are higher order in £/\a\ and V s and 
Vt project onto the 1 Sq and 3 Si partial waves, respectively. Because of Galilean invariance, 
the interaction can only depend on the relative momenta of the incoming and outgoing 
particles k and k'. In a momentum cutoff scheme, the potential in Eq.Q can be regularized 
by multiplying with a regulator function, exp[— (k 2 + k' 2 ) n /A 2n ], where n > 1 is an integer. 
The unphysical cutoff parameter A is arbitrary and strongly suppresses the contribution 
of momentum states with k, k! > A. Expanding the regulator function shows that the 



coefficients A 2 \ n of the momentum expansion of the effective theory at order (k 2 + k 
are modified by the regulator. One can either absorb these additional contributions in 



/2\n 



the definition of the couplings A 2 '! 



nly 

2 



terms at orders 
In the end, all 



, 2n or choose a value of n such that onl 
higher than the desired accuracy in the momentum expansion appear 
observables will be independent of A up to higher order corrections. At leading order, we 
can choose a Gaussian regulator with n = 1 and the regularized potential takes the form 



(k'lVbolk) = V s \ s 2 g(k')g(k) + V t X\ g(k')g(k) 



(2) 



with g(k) = exp(-k 2 /A 2 ). 

The interactions in Eq. (J2J) are separable and thus, the two-body problem for each par- 
tial wave can be solved analytically. The two-body t-matrix can be written as t t<s (E) = 
\g)T ttS (E)(g\, where E denotes the energy. The two-body propagator T t)S (E) is given by: 



T t ,s(E) = 1/A*' s -4tt / dqq 



i 9{qf 
E-q 1 



(3) 



Here and in the following we use units with h = m = 1 for convenience. The coupling 
constants A 2 and A 2 can be fixed by demanding that the triplet and singlet scattering lengths 
a t and and a s are reproduced correctly by the corresponding t-matrices. 
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The properties of the triton are determined by the Faddeev equations 23]. To leading 
order in £/\a\, all internal orbital angular momenta can be set to zero. Thus, the triton 
spin 1/2 is built up from the spins of the nucleon only. The triton wave function can be 
decomposed into Faddeev components. We can eliminate all but one of the components and 
obtain an equation for the remaining component ipi : 

fa = G tPfa + G tG t 3 (l + P)fa , (4) 

where G is the free three-particle propagator and P = P13P23 + P12P23 is a permutation 
operator that generates the omitted Faddeev components; Py exchanges particles i and j. 
The full wave function \l/ can be recovered in the end from the equation 

* = (1 + P)fa + G Q t 3 (l + P)fa. (5) 

The auxilliary quantity £3 is the solution of a Lippmann-Schwinger equation with a SU(4)- 
symmetric three-body contact interaction 

(ui U2IV3I111 u' 2 ) = V a X 3 h(u 1 ,u 2 )h(u' 1 ,u 2 ) (6) 

only. Here V a denotes the projector on the total antisymmetric three-body state with total 
spin 5 = 1/2 and total isospin T = 1 /2 as for example given in |24[ . The regulator function 
/i(ui,U2) = exp(— (u\ + §w|)/A 2 ) is defined in terms of the familiar Jacobi momenta of 
the three-body system. The three-body coupling A3 must be determined from a three- 
body observable. This interaction is required in order to renormalize the three-body system 
and achieve independence of the cutoff parameter A. Its renormalization group behavior is 
governed by a limit cycle k| 3| • For large values of A the running of the coupling constant 
is described by 

c_ sin(g ln(A/L 3 ) - arctan(l/s )) 

31 ~~ A 4 sin(s ln(A/L 3 ) +arctan(l/s )) ' 

where c ~ 0.016 is a normalization constant and So ~ 1.00624 is a transcendental number 
that determines the period of the limit cycle. If the cutoff A is multiplied by a factor 
exp(ri7r/s ) ~ (22.7) n with n an integer, the three-body coupling A 3 is unchanged. The 
three-body parameter L3 can be determined directly from observable quantities like the 
triton binding energy B t . 

Since all phenomenological nucleon-nucleon potentials are fitted to the same two-body 
data, they differ only in their off-shell properties. It is well known, however, that short- 
distance three-body interactions and off-shell two-body interactions are related via unitary 
transformations (See, e.g., Ref. j25| and references therein). As a consequence, the difference 
between the various phenomenological nucleon-nucleon potentials corresponds at leading 
order in £/\a\ to different values of the three-body parameter L 3 . The variation of L 3 
parametrizes the Phillips and Tjon lines (l3l lil 16| . 



Having expressed all relevant coupling constants in the effective interaction potential in 
terms of physical observables, we are now ready to compute the triton charge form factor. 

III. TRITON CHARGE FORM FACTOR 

At leading order (LO) in £/\a\, it is sufficient to apply the electric charge operator to 
the wave function in the impulse approximation. At next-to-leading order (NLO) the linear 
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effective range correction has to be included as well. For neutron-deuteron scattering, various 
higher order calculations have been performed |l2|, HfJ 27, |28[ 2^|, but no calculations 



are available for the triton charge form factor. We note that using minimal substitution 
in the momentum dependent potential (J2J) will generate additional contributions from the 
regulator functions which raises the issue of gauge invariance. However, as in the case of 
the momentum expansion of the effective potential discussed above these terms are at least 
of order k 2 /A 2 . They can either be absorbed in counterterms of the corresponding order 
(if present) or shifted to higher orders by choosing regulator functions with sufficiently high 
powers of n [22]. As a consequence, the effective theory calculation will be gauge invariant up 
to higher order corrections suppressed by (k 2 /A 2 ) n . In this paper, we focus on the universal 
aspects of the triton charge form factor and work at leading order. Gauge dependent terms 
would enter at order k 2 /A 2 which is beyond the accuracy of our calculation. Higher order 
corrections will be considered explicitly in future work |30j]. 

A part of the following calculation has to be performed in configuration space, thus 
we define the corresponding Jacobi coordinates x = r x — r 2 , y = r 3 — |(ri + r 2 ), and 
R = |(ri + r 2 + r^). The charge form factor of the triton is given by 

F c (q 2 ) = (^K+qk/lpclki^K) , (8) 

where q = kj — kj, kj and ky are the initial and final momentum of the scattered electron, 
and \&k denotes the full wave function with center of mass momentum K. The charge 
density operator pc is defined as j3l| 

Pc 



+ r iz )ff c (r - r,) + \{l - r iz ) Pc {v - Ti )] , (9) 



where p c and fP c denote the charge densities of the neutron and proton in configuration space, 
respectively, and are related to the familiar Sachs form factors by Fourier transformation. 
The rj denote the positions of the constituent nucleons. At LO the structure of the nucleons 
does not contribute and the charge form factors of the proton and neutron are simply G P E = 1 
and G% = 0. The first correction from nucleon structure comes from the nucleon radii and 
would contribute at N 2 LO in our counting scheme. 

Due to the symmetry of the wave function, we can set pc = 3^3. Furthermore, we can 
drop constant overall factors as we will normalize Fq at q 2 = to unity in the end. Thus, 
by inserting a complete set of states into Eq. (JHJ) and transforming to Jacobi coordinates x, 
y, and R, we obtain 

F c (q 2 ) = /d 3 xd 3 yd 3 i?^(^K +q |xyRi)p 3 (q)exp(z^qy) 

J i=l,2 6 

x exp(iqR)(xy Ri|^ K ) • (10) 

Since the full wave function factorizes, (^K|xy R) = (^jxy) ■ exp(— iKR), we can elim- 
inate the center-of-mass integration. Fourier transforming to momentum space, we finally 
obtain 



;q 2 ) = / dV d\ 2 ^(^| Ul u 2 i)p3(q)( Ul ,u 2 - -q, . (11) 
J i=l,2 6 



where ui and u 2 are the respective conjugate momenta. Using Eq. (jSJ), the form factor can 
be expressed in terms of integrals over the first Faddeev component ipi. The corresponding 
expressions are somewhat lengthy and can be found in the Appendix. 
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FIG. 1: The triton charge form factor. The solid (dashed) line denotes the leading-order result 
using at and a s (Bd and aj as two-body input. The black circles and squares are experimental 
data from Refs. [32( and |33(], respectively. 



IV. RESULTS AND DISCUSSION 



We are now in the position to calculate the charge form factor of the triton. A simple 
estimate of higher order corrections can be obtained by comparing the results from using 
the experimental values for the triplet scattering length at and the deuteron binding energy 
Bd as two-body input. To leading order in £/\a t \, these two quantities are simply related 
via Bd = 1/af; at higher orders they differ by effective range terms. In Fig. Q we show the 
results for the triton charge form factor at LO. The solid line denotes the calculation using 
a t and a s as two-body input, while the dashed line denotes the calculation using B d and a s 
as two-body input. In both cases the three-body parameter L 3 is chosen to reproduce the 
experimental value of the triton binding energy B t = 8.48 MeV. The corresponding values 
are L3 = 21.2 fm _1 for the solid line and L3 = 19.8 fm _1 for the dashed line. The form 
factor results are plotted up to q 2 = 0.3 fm -2 . For momentum transfers of the order of 
the pion mass, the long-range character of one-pion exchange becomes important and the 
pionless theory breaks down. 

The black circles and squares are experimental data from Refs. |32| and [33| , respectively. 
While the general trend of the low-momentum transfer form factor data is reproduced by 
our calculation, we obtain a somewhat larger slope than the experiments. 

The slope of the charge form factor at low momentum transfer defines the triton charge 
radius: 

F c (q 2 ) = l-q 2 (r 2 >/6 + ... . (12) 

For notational simplicity, we use rc = (r 2 ) 1 / 2 in the following. 

Fitting a polynomial in q 2 to our result for the form factor, we can extract the triton 
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charge radius. We have fitted polynomials of varying degree up to 5th order in q 2 and 
verified that the extracted radius is independent of the degree of the polynomial. This 
procedure leads to a triton charge radius rc = 2.05 fm for a t and a s as two-body input and 
r c = 2.18 fm for B d and a s as two-body input. The expected error of our LO calculation is 
£/\a\ ~ 30%. Averaging the results from the two equivalent choices for the two-body input, 
we obtain r c ° = (2.1±0.6)fm. Our result is about 0.4 fm larger than the experimental value 



r c XP = (1-755 ±0.086) fm |3J] but both values are consistent within their error bars. It would 
be valuable to extend the calculation to higher orders in order to judge the convergence of the 
expansion in £/\a\ |3(|. This requires including the linear effective range correction at NLO. 
Higher order one-body terms from the nucleon radii contribute only at N 2 LO. Inclusion of 
the nucleon radii from Ref. [35| would increase our result for rc by about 0.1 fm. The 
magnitude of this correction is consistent with our error estimate of 0.6 fm at LO. We note, 
however, that there are additional contributions at N 2 LO such as a counter term that arises 



from gauging the momentum-dependent three-body force found in Ref. |2j 

In Ref. |l9j, Friar, Gibson, Chen, and Payne observed a correlation between the triton 
charge radius and binding energy. They explained this correlation as resulting from the 
primary sensitivity of radii to the outer parts of the wave function which are determined 
by the triton binding energy. Using the asymptotic form of the S-state wave function in 
hypersperical coordinates, ip oc exp(— p\^Bl)/ p 5 ^ 2 , they were able to predict the leading 
functional dependence 

r c = bB t - 1 ' 2 . (13) 

If the asymptotic form of the wave function is used for all hyperradii p, one obtains 6=1/2 
[T9| . The correlation curve from Ref. 0] is reproduced by the phenomenological expression 

/ D \ -0.41 



In Fig. El we show the correlation between the triton charge radius and binding energy at 
leading order in £/\a\. The solid (dashed) line denotes our calculation using a t and a s (B d and 
a s ) as input parameters. Both lines are parametrized by the three-body parameter L 3 . The 
circles indicate Faddeev calculations using different potentials from Ref. [3^, |3^ while the 
correlation curve of Ref. H| is given by the dotted line. The square gives the experimental 



values. The two triangles show modern calculations (including meson exchange currents) 
based on the AV18 potential with and without the Urbana IX three-body force j38|. The 
potential models give the same two-body physics but the corresponding values of L 3 are 
generally different. As a consequence, the observed correlation between the triton charge 
radius and binding energy is naturally explained by variation of L 3 . It has the same origin as 
the Phillips and Tjon lines. We note that the dependence of the calculated correlation curves 
in Fig. |21 on the triton binding energy is similar to Eq. (|13|) while the prefactor is different. 
Using the same units as in Eq. (|14|) . we find a prefactor of 7.6 (9.5) and an exponent of 
—0.61 (—0.69) for the solid (dashed) curves. This leads to a slightly stronger B t dependence 
in the leading order effective theory as in the model calculations. 

In the case of the Phillips line, the inclusion of the range correction moves the universal 
line closer to the experimental values 0, |27}. Since the linear effective range correction 
is the only new contribution at NLO, we expect a similar improvement for the correlation 
between rc and B t . A definite conclusion, however, requires the explicit calculation of the 
range correction which is beyond the scope of this study. 
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FIG. 2: The correlation between the triton charge radius and binding energy. The solid (dashed) 
line denotes our leading-order result using at and a s (B^ and a s ) as input parameters. The circles 



indicate Faddeev calculations using different potentials from Refs. |36l . 1371 ] while the two triangles 
show modern calculations using AV18 with and without the Urbana IX three-body force ji^]. The 
correlation curve of Ref. 19] is given by the dotted line. The square indicates the experimental 
values. 



It is also interesting to note that there is no universal relation between the binding energy 
and the position of the minima and maxima of the form factor or the magnetic moments 
[ioj . For the latter quantities, the universality is destroyed by meson exchange currents. In 
the language of the pionless effective theory this corresponds to two-body currents entering 
at low orders in the expansion in £/\a\. It would be very interesting to study the mechanism 
that destroys universality for the magnetic moments. The position of the minima and 
maxima of the form factor (q 2 > 14 fm~ 2 ), however, is outside the range of validity of the 
pionless effective theory. 
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APPENDIX A: FADDEEV COMPONENTS 



In the three-body system, it is always possible to choose a cutoff A for which the three- 
body force vanishes. Thus, the three-body force contribution can be dropped in the following. 
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Using Eq. (jSJ), we express the full wave function \1/ in terms of the first Faddeev component 
tpi and work in a partial wave projected basis. Defining q = \u 2 — fq|, we have 



du 1 u 2 1 du 2 ul j ds|(^|n 1 u 2 l)p 3 (g)(n 1 gl|(l + P)|^i) 

+{^\u 1 u 2 2)p 3 (q){u 1 q2\(l + P)\tf; 1 )} . (Al) 
Applying the usual overlap equalities we obtain 

^c(q 2 ) = l duiu\du 2 u 2 2 I dx\ (^i|«iw 2 l)p 3 (g)(wigl|T/>i) + (i/) 1 \u 1 u 2 2)p 3 (q)(u 1 q2\'4> 1 ) 



dx' 



-i 



1 .1 

-(V>i|«i(itl, u2) u 2 l)p 3 (q)(uiq l|-0i) + -(ipi\ui(ul, u2) u 2 2) p 3 (q) (mq 2|-0i) 



+^(0i|ui«2 l)p 3 (g)(Mi(wl, ?) ^2 + ^(^i|«iM 2 2)p 3 (g)(ui(nl,g) m 2 2|-0i) 



i 

-i / dx'drc" 



tt2) m 2 l)p 3 (g) («i(itl, g) u 2 



+ w2) u 2 2)p 3 (q)(u 1 (ul, q) u 2 2|-0i) 



Now we separate the spin- and isospin components from the momentum space function 



|1) = X1V2 , 
1 2) = X2V1 , 

where \ denotes the spin- and rj the isospin part and 

(wiw 2 l|^i) = (mi«2|V ;1 )Xi'72 
(u x u 2 2\il)i) = (iHU^ip^XiVi 

The isospin functions are written as 



1 
1 

7! 



+ + - + 



2(+--) 



Thus, if we apply p 3 on the isospin functions we get 

p 3 (q) Vl = -L[G p E (q)(--+)+Gl(q)(- + -)-2G n E (q)(+- 



Ps(q)V2 = -j=\G p E (q)(- 



G E (q)(- 



(A2) 
(A3) 



(A4) 

(A5) 
(A6) 

(A7) 
(A8) 
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where G P E and G E denote the proton and neutron charge form factors, respectively. Inserting 
the above definitions into Eq. (jA2|) we obtain 



G P E (q) [±l}(q) + ll?(q) + g#(g) + ^if(g) 

+lim+^(q)+^ 6 m+^) 

+G n E (q) [\il(q) + \iKq) + \il{q) + ^"(q) 



(A9) 



where we have used the definitions 



I\ 



dui u\du 2 ul / dx \u\Ui) {uiq\ip\) 



du 1 ujdu 2 ul I dx J dx'(i)\\u 1 (u 1 ,u 2 )u2(u 1 ,u 2 ))(u 1 q\'ip 1 1 ) 
duiu\du 2 u\ I dx I dx' (ip\\uiu 2 ) (iii(ui, q)u 2 (ui, q)\ipl) , 



-i 



-i 



duiuldu 2 ul / da; / da;' / da;' / dx" (ipWui^i, u 2 )u 2 (ui, u 2 )) 



-i 



-i 



x(tti(ui,g)«2(ui,g)|Vi) • (A10) 

At leading order, only the total charges of the proton and neutron contribute in Eq. ()A9j) : 

G p E (q 2 = 0) = 1 and G™ (q 2 = 0) = . (All) 

The first correction from nucleon structure is due to the charge radii of proton and neutron 
and enters at N 2 LO. 
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